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The DeWitt expansion of the matrix element 
M xy = (x\ exp — \\{p — A) 2 + V]t \y), (p = —id) in powers of t can be made 



o 

On 



in 

On 



in a number of ways. For x = y (the case of interest when doing one-loop 
calculations) numerous approaches have been employed to determine this ex- 
pansion to very high order; when x ^ y (relevant for doing calculations beyond 
one-loop) there appear to be but two examples of performing the DeWitt ex- 
Qh! pansion. In this paper we compute the off-diagonal elements of the DeWitt 

expansion coefficients using the Fock-Schwinger gauge. Our technique is based 
on representing M xy by a quantum mechanical path integral. We also gen- 
eralize our method to the case of curved space, allowing us to determine the 
DeWitt expansion of M xy = (x\ exp ±[-^(d M - iA^gT \fg{d v - %AJ)]t\y) by 
use of normal coordinates. By comparison with results for the DeWitt expan- 
sion of this matrix element obtained by the iterative solution of the diffusion 
equation, the relative merit of different approaches to the representation of 
M xy as a quantum mechanical path integral can be assessed. Furthermore, 
the exact dependence of M xy on some geometric scalars can be determined. 
In two appendices, we discuss boundary effects in the one-dimensional quan- 
tum mechanical path integral, and the curved space generalization of the 
Fock-Schwinger gauge. 
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I. INTRODUCTION 



There is a long history of computing the elements a n (xo, A) in the expansion of M xy in 
powers of t, 

M xy = (x\ exp {- [\{ V - Af + V] t} \y) (la) 

-A 2 /2t oo 

= (2^S OB(iB0 ' A) * B ' (lb) 
(£> = no. of dimensions, p = — z<9, x = (x + f/)/2, A = x — y) , 

in the limit A = [|TJ. This expansion when A = is extremely useful when examining 
certain properties of the generating functional at one-loop order; in particular, the divergence 
structure of a theory at one-loop order can be discerned. Among the approaches used 
to evaluate a ra (x ,0) are the perturbative solution to the heat equation 0^], the use of 
pseudo-differential operators J3|, working in momentum space |3J, systematically rearranging 
a Schwinger expansion of ([la]) in powers of A and V into an expression of the form (|lb|) 
|5||| and representing flTa]) as a quantum mechanical path integral (QMPI) hence expanding 



it in powers of t |7||§. The only places of which the authors are aware where a n (xo, A) is 
considered for A ^ are in M and ||. These coefficients are useful in considering multi-loop 
processes [[J, which motivates us to pursue them further. The quantum mechanical path 
integral has proved useful in computing Green's functions at one-loop order [|TT]|-|T3|] and 
beyond |T^-[16|; this suggests using this approach to examine a n (x , A) for A^O. Although 
our method is not identical to that of [^], the two approaches are similar and both results 
agree when A = 0. 

The representation of 

M xy = (x\ exp [\ [j={d, t - %A^^fg{d v - iA v )\ t} \y) 
in terms of a quantum mechanical path integral is not uniquely specified []T7Hl9[1 , as discussed 



in 



T4| . We use one of the various forms of the QMPI to expand M xy and compare our results 
with those of ||. Furthermore, a partial summation of the DeWitt expansion to obtain the 
full dependence of M xy on R and -R^A^A^ is possible [ pOf . 



II. EXPANDING M X Y 

It is possible to represent M xy as an appropriately normalized QMPI [f2T 



M xy = J*Dq(T)Vex V jT* dr j-^M + iq{r) ■ A(q(r)) - V(g(r))| (2) 

where path-ordered integration is implied over trajectories with end points g(0) = y and 
q(t) = x. 

We attempt to construct a power series about some point xq which we arbitrarily choose 
to be the mid-point between x and y. Defining the relative coordinate 5 by 
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q(r) = x + 8(t) , 
and imposing the Fock-Schwinger gauge condition [f2!| , 

8(t)-A(x + 6(t)) = 0, 
one can expand the gauge field in powers of 6, 

A^xq + 8(t)) = [ daa5 x (r)F Xf ,(x + aS(r)) 



OO 



[5(r) ■ D(x )f 8\r)F x ,(x ) 



The scalar potential can be similarly expanded, 



V(x + S(t)) 



Af=0 ly ■ 



D(x )] N V(x ) 



Here gauge- covariant differentiation at x has been denoted by D(x ). Together, (|5 
(|) allow (|la|) to be written as 



M. 



XIJ 



A/2 
-A/2 



D5{t) exp 



o 



OO 1 

J L=0 n - 



(3) 
(4) 

(5a) 
(5b) 

(6) 
and 

(7) 



I 00 1 



dr[5(r) ■ D{x )\ 



N 



N + 2 



5»(t)5\t)F x ^x ) - V(x ) 



In the above expansion we intend to treat all terms in the expansions of the potentials as 
perturbations on the free-field action | Jq <ir5 2 (r). By contrast, the authors of reference |23| 
have shown how, through the introduction of an appropriate tensor basis, one can include 
the lowest-order term of the derivative expansion of the electromagnetic field in a non- 
perturbative fashion. Although the latter technique could be applied here, we adopt the 
purely perturbative approach which is algebraically simpler and more suitable for the purpose 
of illustration; furthermore, we are free to consider a non-abelian gauge group. 

The path integral in (0) can be evaluated by systematic functional differentiation of the 



standard result [pl|,f24|,p 

rA/2 



/■A/2 rt 

I D5(t) exp / dr 

J-A/2 JO 



5 2 (t) 



+ 7 (t) -5{t) 



(8) 



e -A 2 /2t 

(2nt) D / 2 



exp 



dr 



A • 7 (r) 



dr dr'G(r, t')j(t) ■ 7(r') 



T — T 



Hr + r') + 



TT 
t 



with respect to 7 a (r) and then setting 7 = 0. (Here, G(r, r') = 

is the Green's function of a free particle on the worldline.Q) For example, after two such 
derivatives, it is easily shown that 



1 In appendix [A| we address some concerns about the validity of this Green's function on the finite 
interval [0, £]. The complications presented there are not expected to contribute in this flat-space 
limit. 
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A/2 



DS(r)8 a (r a )8^T b ) 



-A/2 



exp 



Jo 2 



(9) 



T, 



t 



1 

2 



From (P), it is easily seen that no term in (^) will involve factors of A 2 . From this 
observation, combined with simple power counting arguments, it is straightforward to tab- 
ulate the possible contributions to the various coefficients a„(xo, A) (see table |). (When 
a temperature dependent QMPI is considered as in |12[ , then the temperature provides a 
second dimensionful parameter which must be considered). The coefficients of the various 
contributions can be easily determined by appropriately choosing L and N in (|7]) and then 
systematically applying (|]). For example, if L = n, N = 0, then it is apparent that the 
contribution to a n (xo, A) proportional to V n (xo) is -^ (— V(xo)) n - By setting L = N = 1 in 
(0), we find after a very short calculation that the contribution to di(x , A) proportional to 
A a DpFj is ^A a D p F a t3 . With L = 2, iV = 0, the contribution to M xy is straightforwardly 
computed to be 



e -A 2 /2t 
(27Tt) D / 2 



/ t 2 N 

2^ r an r p 48 



giving a contribution to both ai(xo, A) and 02(^0, A). These results are all consistent with 
the flat space limit of the expressions for a n (x , A) given in 0. 



III. EXPANDING M X y 



As has been noted in the introduction and in fl4| , there are various representations of 
the matrix element 



M, 



■1:1/ 



x\ exp {I [-Wfy - iA^Jgidv - iA v )] t) 



(10) 



in terms of a QMPI. We adopt the approach of [|T^] in which the QMPI is computed using 
a normal coordinate expansion of the coordinate being integrated over and the form of 
the classical action is fixed by having M xy satisfy the appropriate heat kernel equation. 
(This representation does not coincide with the expression given in |25] for any value of the 



parameter p appearing there.) This representation gives a dependence of M xy on R that 
coincides with that of [^,0, and agrees, to all orders so far checked, with the results of 
The full dependence of M xy on i^^A^A^ can also be determined. 
We are now faced with evaluating 



M. 



Dq a (T)Jg(q(r))Vexp 



(11) 



The factor J g(q{r)) = det 1 ^ 2 g^(q(r)) in the measure also occurs in the non-linear sigma 



model PB], but there it is usually discarded as it gives a contribution to the effective action 
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that is proportional to 6(0) which, when regulated using dimensional regularization [27 
or operator regularization ||28|| , goes to zero. We are dealing with a model in which no 
regularization is required; indeed it turns out that divergent contributions from ^fg are 
essential to render the path integral in (|ll|) well-defined. It is most convenient to incorporate 
the effects of y/g by using ghosts as was done in . Perhaps the simplest way to do this is 
to introduce a vector of real Bosonic fields b a (r) which vanish at the end points (r = 0, t). 
(The structure of the ghost sector is not unique, and the authors of |Ll| opt, instead, to use 
either a pair of Fermionic scalars or a contravariant Bosonic ghost b a and two contravariant 
Fermionic ghosts c a and c a .) Then, equation ( |TT] ) can be re-expressed as 



xy 



jT Dq a (r) f ° Db p {r)V exp { - f\h 



\g^b,(r)K(r) 



o l Jo 



A normal coordinate expansion PS| is now made about a point 0(r) so that 

Following the developments of |27|, we have 

R(q(r)) = i?(0(r)) + i£ Q (^(r))f (r) + ^R, a ^(r))^ (t)^ (r) + • • • , 



(12) 



(13) 



(14a) 



9,u(q(r)) = 9^(r)) + \R^{4>{r))i° '(r)^ '(r) + \R m p^(r))C (r)^(r)C(r) (14b) 
+ (iiWry*(0(T)) + 1 5 R^Mt))R%M{t))) C(r)e(r)C(r)e(T) 
+ ■■• , 



gTW)) = <T{m) - |^V(0(r))r(r)^(r) - |^/ ;7 (0(r))r(r)^(r)e(r) (14c) 
+ (-^V^^(r)) + J^% ct (0(t))#V(^))) r(r)^(r)r(r)^(r) 



+ 



and 



g"(r) = r{r) + £> t £"(t) + |i?W0(r))r (r)^(r)0 7 (r) + 



(14d) 



As shown in appendix [FJ, by imposing a gauge condition analogous to (^), one finds that 
the corresponding normal coordinate expansion for the gauge field is ( p3|) 



+ (jiWr(^W) + ^^-(0(r))^^(0(r))) ^WfW 
+ (£iWy*(0OO) + ^F aCT (0(r))i?^ w;5 (0(r)) 

+^^ ;5 (0(r)) J Rv,(0(^))) rwevjrwev) + 



(14e) 
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If we take 4>(t) to be the geodesic mid-point, x$, between x and y, then 4> a (r) vanishes and 
the above expansions simplify a bit. Letting A denote the difference between the normal 
coordinates of x and y, equation (|T2| ) becomes 



M 



A/2 f 

DC(r) / Z%(r)exp 

-A/2 JO 



dr \y^e(r)t(r) + \g^b,{r)b u {r) 



-Rt/8 



(15) 



00 1 f ft 

AT=0 iV • 1 J ° 



-i (r) + lF aKP e(r)^(r) + ■ • •) £ 



WW 



A" 



(All geometrical and gauge quantities in (p~5| ) are evaluated at x .) 
The standard results 

A/2 D^explfdr (-\g^{r)i v {r) + ^{r) lii {r j 

-A/2 UO v 

exp \l n dT (-- + -) A" 7 „(t) - § / n rfrrfr'G(r, tV7»7v« 



(16a) 



and, 



J° Db a (r)ex P y*dT (-\ g ^b,{r)b v {r) + & m (t)B"(t))} 
= v^expji J* dTdT'G^ ost (r,T') gtlu B^r)B^r' 



(16b) 



B(0)=B(t)=0 



permit one to compute the functional integrals appearing in (15). Equivalently, one can 
generate the necessary contractions using the formalism of equation (|A2j). In either case, 
one should be careful to use the Green's functions of equations (|A6|) (or in simple cases, 
equations (|A8|) ) as discussed in appendix |A|. 

For example, if we restrict our attention to the contribution to (|1|) that is linear in 
R^avp, we have 



M. 



R 

xy 



A/2 
-A/2 



DC(r) / Db p {r)e 



-Rt/8 



x exp 



dr 



(2nt) D /< 



e -Jtt/8 /" rfr j _ 1 r^a/3 G ( f; ^g,^ r ) + { g m g vP + gtfgVajGtfa f ) 



-A M A 



"(|)Vg(t, 



r 



H + l) G(f,f) 
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- (A>*A a g ul3 + A"AV tt + A^AV" + A u A a g^) \ (-\ + f) G(t, t) 
-A a A p A"A v (-§ + £ 



+ 



-G{T,T)sf* + (-\ + z) A a A? G^ host (r,r)g^\[-R 



(17) 



Here, a dot over an argument of the Green's function G indicates differentiation with respect 
to that argument. As explained in appendix 0, these derivatives cannot always be obtained 
at the end points (r = 0, t) by naive differentiation of G(r,r') (c.f. equation ( |A8cj )). When 
the explicit form of G(t,t') and its derivatives are substituted into we are left with 
with 



R 

xy 



e - gtlv A»A"/2t i 

(2nt) D / 2 6 C 
R a ?A a AP 

+R 

e - 9tl „A»A»/2t 

(2nt) D / 2 



-Rt/8 



dr 



.1 + 1 

2 ~ t 





'-T H 




-H- 













- [6{2t) + 5{2t - 2r)] 
[5{2t) + <5(2t - 2r)] (-t - 



2 ~ t 



-Jtt/8 



(18) 



where i? Q/ g = R^aPn and i? = R a0 g a ^ . (Notice that all dependence on 5(0) in fllTD has 
cancelled out due to the compensating contributions from the ghost-fields.) It is easily seen 
that the entire dependence of M xy on R and R a pA a A 13 is given by 



jCfiRA-R-A) 

xy 



e -g^A^A"/2t oo e -Rt/8 



(27Tt) D / 2 
- gilv A^A»/2t 

{2llt) D l 2 



E 

iV=0 

exp 



R a0 A a A p Rt 



N\ { 12 
_m _ R a0 A a AP 
~12 12 



+ 



24 



TV 



(19) 



The dependence of M xy on both R and i? Q , / gA Q A /3 in ( ]T9| ) agrees with that of p0|j2]fl9|1 (once 
the different normalization of t is taken into account). 

A completely analogous calculation can be used to fix the lowest order dependence on 
.Rc^A^A^A 7 ; it is found to vanish. This also appears to be consistent with the results of 
(where the coefficient -R Q/ g ;7 is evaluated at the end-point y, instead of the mid-point Xq). 
Further terms in the DeWitt expansion of M xy can be similarly determined. 

The techniques used in this section may be employed to find the effective action for 
a particle moving in a gravitational field [ SO] . This involves taking </> a (r) to be arbitrary 
in (|T3|) rather than restricting it to be Xq. 



IV. DISCUSSION 

In the preceding sections we have considered how the QMPI can be used to determine 
the elements of the DeWitt expansion for the heat kernel both on and off the diagonal. This 
technique is seen to be easier to use than the original approach in which the heat equation 



7 



was solved perturbatively. By employing the off-diagonal elements, calculations can be done 
to two-loop order ||. The method works in both flat and curved space. 

Because of the simplicity of the proper-time (r) integrands which arise when implement- 
ing this method, we expect that high-order covariant expansions could be computerized 
using presently available symbolic algebra packages. This approach has already been em- 
ployed in the diagonal case j7|. In the off-diagonal case, surface terms (proportional to A) 
introduce an additional combinatorical consideration which should be tractable. 

Finally, it is worth noting, that although our treatment of this problem focused specif- 
ically on the case where the vector potential was a gauge field, this is unnecessarily 
restrictive and the methods presented here would work equally well for a general vector cou- 
pling to a background field. For example, if one were interested in radiative corrections to 
Fermionic Green's functions in flat-space QED, the following quantum mechanical operator 
would be of interest |yj 



p-'.i^v - (1 - J) PuPv -(i^) T {ihn) ^ 

(7i/0) j> — m 

V -($lv) T -{j>- m) T j 



(20) 



By factoring in the constant supermatrix 



QvX 
(j> + m) T 

-ty + m) 



(i/,A) 



and choosing the a = 1 gauge, operator (|20| ) becomes 



/ P 2 9ti\ -07 M )(^ + m ) -(ri^) T {.i> + m) T 

{l\4>) p 2 + m 2 

v. -(^7a) T 



(21) 



o 9 

p + 771 



(M,A) 



After completing the square of p (and noting that p T = —p in the coordinate-space repre- 
sentation), the heat kernel of this operator is easily shown to have the form of equation ( la ) 
with vector potential, 



AM) 



( ^(qh^lu -i.lvl^{q)f 

















(22) 



and scalar potential 
/ 

V(q) = 



MM 2" 
2 



(^(g) 7f 
T 



— 771 



o 



V -(V»(?)7a) 



m 



) 



After the appropriate Taylor expansions of these potentials is substituted into eq. 
method should proceed in the obvious way. 



(23) 



the 
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APPENDIX A: EVALUATION OF THE PATH INTEGRAL 



Here we explicitly evaluate the path integrals used in the previous sections. Firstly, we 
define the type of integrals which need to be evaluated to obtain the expectation value of a 
function F of dynamical variables £ m (t) and 6^(r), 



(F(^6)) A = (f,0 F(£,b) 



-± 

2 ' w 



{27lt) D / 2 



A/2 
-A/2 



DC(r) f Db p {r)F{i,b) 
Jo 



(Al) 



x exp 



dr 



which satisfies the normalization condition (I) a — 1- 

We eliminate the boundary parameter A from the path integral by integrating over 
fluctuations about the classical geodesic £ci( r ) = ( — 1/2 + r/t)A. Letting £ — > £ + £ c i it is 
easily shown that 



(F[{,b)) A = (F[Z + £ ch b)) . 



(A2) 



Thus, for the purpose of evaluating any particular term in equation (|T|), it will be sufficient 
for us to concentrate on evaluating (F[£, b])o where F is a monomial in £, £ and b. We can 
follow the procedure usually used in applying the path integral in field theory and so we 
only need to evaluate the various two-point functions, (6 At (r)6 I/ (r / )), (^(^^(r')), etc, and 
then apply the appropriate Wick expansions. 

Jumping directly into the continuum limit, one would obtain, using the standard tech- 
niques pp 



(e(r)C(r')) = -G(r,r')g 



(A3a) 



T 



t'\-Ut + t') + 



TT 



(b,(r)K(T')) = G^ st (T,T')g^ 
= 6(t 



(A3b) 
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{e(r)K(r')) = 



(A3c) 



Upon close inspection of the above propagators one finds that the method used for finding 
the explicit form of the Green's functions is problematic, especially at the end points. For 
example, the defining equation for G 



Or 2 



G(r,r') = 5(t-t') 



and the corresponding boundary conditions 

G(t, 0) = G(r,t) = 0, 

are contradictory at the r' = and r' = t boundaries. Furthermore, the explicit form of the 
ghost propagator of equation ( |A3b|) does not satisfy its homogeneous boundary condition 
at the corners (0,0) and (t,t). 

In order to circumvent these difficulties, we adopt the approach of [IS] where the path 
integral is taken to be the limiting case of integration over a finite set of discrete Fourier 
modes. The Fourier expansions 



M 



TVKT 



M 



sin 



n=l 



t 



TVKT 



sm 



71=1 



t 



(A4) 



automatically incorporate the necessary boundary conditions in the A = path integral; 
furthermore, they ensure that the path integral is over functions which are periodic, as 
required in [10]. If the cutoff M is finite then the path integral over Fourier coefficients 
is well defined. Substituting these expansions into the action of (A~I), we find that the 
propagators of the Fourier modes are simply |L9 



\SnSm/ 



2t 
n 2 7T 2 



3nm9 



^ unburn) . ^nm9fj.u 



(A5) 



which we take to be our fundamental propagators. Then using (|A4|) the field propagators 
are found to be 



M 



(e(r)c(r')} = grtY, 

n=i 
I M 

(6 M (r)^(r)) = g^u- E 
(e(r)K(r')} = 0. 



1 



2^2 



rm{r — r') n7r(r + r') 

cos cos - 



t 



cos ■ 



nulr + t') 



cos • 



t 



(A6a) 

(A6b) 
(A6c) 



These are the propagators which should be substituted into any Wick expansion of (F[£, b]) . 
In principal we should attempt to approximate the continuum limit by letting M — > oo 
only after the proper-time integrals have been performed. Unfortunately, these harmonic 
expressions for the Greens functions can be unnecessarily cumbersome; in practice whenever 
the appropriate product of Green's functions is sufficiently well behaved as M — ► oo, we 
will sum them in that limit before substituting them into the Wick expansion. This is most 
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easily done by using the variables t± = r ± t'. Then, on the relevant intervals r + e [0, 2t], 
t_ G [— t, t], the following limits are easily derived from elementary Fourier theory: 



M 

E 



n=l 



n 2 TT 2 



COS 



E — Sin 



TOT 



1 



n=l 

M 

— sin 

M x 

Et cos 

n=l 6 
M x 

Et cos 

n=l 1 





M- 


->oo 


t 






'ii 


M- 


-»oo 


t 






WKT + 


M- 


->oo 


t 








M- 


-»oo 


t 






WKT+ 


M- 


->oo 



T±\ + 



At' 



_ sgnr _ _ _ r _ 



1 

2t 



5(t_ 

5(t + ) + 5(2t - r + ) 



1 

2f 



(A7a) 
(A7b) 
(A7c) 
(A7d) 
(A7e) 



It is thus found that in many cases the following prescription is sufficient to recover those 
propogators which arise in (p^|): 



(e(r)C(r')) 

(e(r)r(r')) 
(^(r)^(r')) 



-G(r,r') 
-G(f,0 



2 I 
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2 ' t 



9 



host i 



r, r 



2Sgn(T - T 
-8(r - t') + \ - 8(2t -t-t')- 6(t + r' 
[5(r - t') - 5{2t -t-t')- 8(t + r')] ff^ 



(A8a) 
(A8b) 
(A8c) 
(A8d) 



(We note here that (|A8c|) can be derived from ( |A8a|) by requiring that G be extended 
periodically outside of the domain < r, r' < t, i.e G{r + t,r' + t) = G(t,t'). This 
periodicity is akin to the definition of the Green's function given in Jl0| on a compact 
periodic surface. A discussion of this point can be found in |T5|] .) 

Of note is the limiting case where r = r' which appears in the calculation of (|1T[) ; there 
the following considerations apply: 

• The term sgn(r — r') in (|A8b|) is odd in r — r' so its Fourier transform necessarily 
vanishes in the r = r' limit. Thus we should take sgn(0) = (c.f. ( |A7bj )) . 

• In the t — t' limit, the Fourier series of 8{r — r') in (|A8d) and ( |A8d|) leads to a 
regulated representation of 5(0) given by t8(0) = M + 1/2, (c.f. ( |A7d|) ). 



Since the Fourier series of the terms 8(2t — 2r) + 5(2t) (which arise in ( |A8c| ) and ( |A8d| )) 
is even about the respective poles (c.f. ( |A7ej) ), any integrated contribution from these 
terms should be halved as a result of the poles coinciding with the end-point of the 
integration region, i.e. 

jT* dr [5(2t - 2r) + 8(2r)\ f(r) = \ [f(t) + f(0)} . 
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APPENDIX B: NORMAL COORDINATE EXPANSION OF THE GAUGE FIELD 



In this appendix we discuss the construction of a gauge-covariant normal coordinate 
expansion for the gauge potential. 

By analogy with the flat-space case discussed briefly in section [Q| and in refs. [IT|,7|) 



the appropriate gauge condition for this expansion is the synchronous gauge [0 (a curved- 



space generalization of the Fock-Schwinger gauge @ ||22|| ) which fits very well in the normal 



coordinate construction. In the basis of the normal coordinate system, the gauge condition 



is 



fVU0 + vr(O) = O. (Bl) 
Either by integrating along the geodesies (which is formally identical to equation (|5i|)) or 



by using differential forms [32] one can show, in the normal coordinate system, that the 
synchronous gauge leads to a gauge-covariant expansion for the vector potential which 
looks exactly like equation (|)b|) with the gauge-covariant normal coordinate derivative 
D a = + [A*;---]- The latter derivative is not covariant under reparametrization of 
the manifold, however using the methods of reference |Z7J it is straightforward to write such 
normal coordinate derivatives at the origin in terms of the corresponding fully- covariant 
derivatives, denoted by indices trailing the semicolon (;). For example, one can show that | 



D Px F Poi— F Par,Pi (B2a) 

D Pi D Pl F Po7= F Por,PiP2 + I^Ai^ifta (B2b) 

D Pi D P2 D Pi F Poi= F PorAP2P3 + f-^o^xftnifts + F PoS;Pi R p 2 p 3l ( B2 c) 

D Pi D Pi D fi2 D Pl F P01 = F P0TAP2PiP4. + 5 F PoS RS f3lP2rA>P4 + ^O^l-Rft^sTiAl 

+2Fp s ;j 3 1 p 2 Rp 3 p 4l + lFp o£ Rp 1026 R s M41 (B2d) 

where = indicates equality at the origin only after symmetrization of the $ indices. Substi- 
tution of equations ( p2|) into eq. (|5b|) yields the fully covariant normal coordinate expansion 
to fifth-order in the normal coordinates, 



A 7 (0 + vr(O) = \ {F^} e + | {Ppor^A^e 1 + | {Wa + kFfM&Mn} f°C 01 ^ 2 

+ 4T6 { F Par,Pilhp3l34 + |- F /3o<5- R A/3 2 7;/33/34 + 2F PoS;PiRp 2 l3 3r ,f3 4 ( B 3) 

+0 (e) . 



The authors suspect that the fourth derivative of a rank-two tensor implied in reference [27] is 
not entirely correct. The corresponding coefficients presented here for the field strength, equa- 



tion ( p2d| ), have b een verified independently. 
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All coefficients in braces {• ■ •} are evaluated at the origin, where the basis vectors for the 
normal coordinate system coincide with those of the original system. Since the potential on 
the left hand side of this equation is not a vector at the origin, its indices must refer to the 
normal coordinate basis. (This is also true of equations (|14b| )- (|14d|) ). The results of ( |B3"D 
agree with those of to order O (£ 3 ). 
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TABLES 



TABLE I. Possible contributions to the various coefficients of the DeWitt expansion. 



Coefficient 



a (x ,A) 



Contributions 



(A • D) k V, (A • D) k (A a DpF a P) , (A • Pf (A Q A^F/ 



ai (x , A) 



D) k (D 2 V) , (A -P) k V (A • D~f V, 
D) k V (A • D) 1 (A. a D p Fj) , (A • D) k (A a D p Fj) (A • D) 1 V, 
pf (a^A^F/) (A • D) 1 V, (A • Pf V (A • D) 1 (A^A^F/) , 

D f F «v ( A . D y Faf)i 

D) k (A°F QM ) (A • D) 1 (DpFfr) , (A • D) k (PpF^) (A • D) 1 (A a F a ^) , 
pf (A a DpFj) (A • D) 1 (WDpFj) , 
D) k h^DpFj) (A • D) 1 [A a APF ail FA , 
pf }A a APF^FA (A • D) 1 (A a DpFj) , 

pf (A Q A^F/) (A • D) 1 (A a A^ F afl F^) 



a 2 (x , A) 



(A 
(A 
(A 
(A 
(A 
(A 
(A 
(A 
(A 



etc. 



(All fields and covariant derivatives are evaluated at x . 

(M = o,i,2,...) 
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